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Abstract In this paper, we propose using curvatures in digital space for 3D object analysis 
and recognition. Since direct adjacency has only six types of digital surface points in local 
configurations, it is easy to determine and classify the discrete curvatures for every point on 
the boundary of a 3D object. Unlike the boundary simplicial decomposition (triangulation) , 
the curvature can take any real value. It sometimes makes difficulties to find a right value for 
threshold. This paper focuses on the global properties of categorizing curvatures for small 
regions. We use both digital Gaussian curvatures and digital mean curvatures to 3D shapes. 
This paper proposes a multi-scale method for 3D object analysis and a vector method for 3D 
similarity classification. We use these methods for face recognition and shape classification. 
We have found that the Gaussian curvatures mainly describe the global features and average 
characteristics such as the five regions of a human face. However, mean curvatures can be 
used to find local features and extreme points such as nose in 3D facial data. 

1 Introduction 

In recent years, there are many applications related to geometric analysis in 3D image processing, 
especially the methods using topological and geometric invariants. [3] p6] This paper presents a 
preliminary study of using digital curvatures for digital surfaces to classify different 3D objects. 

A 3D object can be represented by one or several closed surfaces (2D-manifolds). Curvatures 
that describe the degree of change at a point on the surface have been used for many years in 3D 
image processing j5] [5] [T5]. The typical technology related to curvatures is the following: (1) to 
use triangulation of the surface, (2)to fit the digital image to a continuous surface (using B-spline), 
and (3) to calculate the standard Gaussian and/or Mean Curvatures. This usually refers to the 
areas of computational geometry. 

On the other hand, methods in digital and discrete geometry have been developed rapidly in 
recent years for 3D computer graphics and 3D image processing [14]. For example, recognition 
algorithms related to 3D medical images can be found in 4J in which Brimkov and Klette made 
extensive investigations in boundary tracking. 

In this paper, we propose using curvatures in digital space for 3D object matching, classifica- 
tion, and recognition. Since direct adjacency only has six types of digital surface points in local 
configurations, it is much easier to determine and classify the digital curvatures for every point on 
the boundary surface of a 3D object. Because of the limited number of local configurations, this 
method may only apply to certain applications. 

2 Background 

Given a set of cloud points in 3D, we assume that they are connected. (Such is the case for 3D 
MRI brain images.) Since cubical space with direct adjacency, or (6,26)-connectivity space, has 



the simplest topology in 3D digital spaces, we will use it as the 3D image domain. It is also 
believed to be sufficient for the topological property extraction of digital objects in 3D. In this 
space, two points are said to be adjacent in (6,26)-connectivity space if the Euclidean distance 
between these two points is 1. 

Let M be a closed (orientable) digital surface in the 3D grid space in direct adjacency. We 
know that there are exactly 6-types of digital surface points (5l [6] . 




Figure 1: Six types of digital surfaces points in 3D 

Assume that Mi (M3, M4, Afs, Mg) is the set of digital points with i neighbors. We have the 
following result for a simply connected M [5]: 

|M3|=8 + |A/5|+2|M6|. (1) 
2.1 Gaussian Curvatures and Mean Curvatures 

If Ki denotes the Gaussian curvature of elements in Mi, i — 3,4,5,6. We have (see discussion in 

m) 

Lemma 2.1 (a) K3 = tt/2, (h) = 0, for both types of digital surface points, (c) = — 7r/2, 
and (d) Kq — — vr, for both types of digital surface points. 

We also have a genus formula based on the Gauss-Bonnet Theorem [5] 

g^\ + (IA/5I + 2 • iMel - iMgD/S. (2) 
If Hi denotes the mean curvature of elements in M^, i = 3, 4, 5, 6. We have 

Lemma 2.2 (a) = (b) H4 = for flat neighborhood, and H4 = V2 for a bended neigh- 
borhood, (c) = I . (d) Hq — 0, for both types of digital surface points. 

This lemma can be proved easily based on the formula derived by Meyer et al in 2002 [T7]. 
Goodman-Strauss and Sullivan used the Hq ~ points to construct cubic minimum surfaces j l2| . 
A minimum surface can be defined as a surface whose mean curvature at every point is 0. Here 
we give a brief proof for Lemma 2.2. 

The formula derived in [17] for the mean curvature normal at point x is: 

H{x) ■ n = ^^S^^gAr(^)(cota,; + cot^i)(a; - x.i) (3) 

where N is the set of (discrete) neighbors of x. Ui and Pi are angles in two different triangles that 
are both opposite to line segment xxi] which is shared by the triangles. A is called the voronoi 
region of x. For digital space, the voronoi region is easy to determine. See Fig 2. 

H{x) -11= ^^.(zN{x){x- Xi) (4) 



Figure 2: The voronoi region in digital space 



Note that H{x) • n is a vector. |(a; — Xi)\ = 1 and cota^ = cot f3i =1. in Fig 2. A = j • i 
for M, points. |Sxi G N{x){x - Xi)\ = ^(3) for Mg; e iV(x)(x - a;,)] = for M4 flat, 
and y(2) for Af4 bended; jEa;, S A^(x)(x - = 1 for M5;\T,x^ G A^(a;)(a; - x^)\ = for Mq; 
Therefore, H3 ^ ttVI^) = For bended M4 points (Fig. 1 (c)), i?4 = ttx/I^) = J{2). 

4 Y (3) '*4 

2.2 Digital Principal Curvatures 

The principal curvatures at a point p of a surface, denoted ki and fc2i are the maximum and 
minimum values of the curve curvature on normal planes (intersecting with the surface). The 
relationship among the principal curvature, the Gaussian curvature, and the mean curvature are: 
K = fci • fca and iJ = (fci + fc2)/2. Therefore [T5][17], 

ki= H+ ^Jm - K, k2 = H - \/H^ - K 

It is easy to get, 

y^f"^ = 0.369675. (b) 
- ! + Vi + f = 2.28687, 
= -0? = -1.77245 for both 

From these lemmas, we can find some interesting information about digital principal curvatures. 
k\'^ does not have tt involved in the fomula. However, from Fig. Ic to Fig. la is just like to make 
a 90 degree angle similar to the case from Fig. lb to Fig. Ic. Where the tt came from? Due to 
the fact of digitaldiscrete mean curvatures were obtained from an approximation, we think that 
the value of the mean curvatures for Mi points can be modified as H' = ci ■ H + C2 as long as we 
keep H'^ > H^^ > > = Q, where ci and C2 are two constants. 

3 Digital Curvatures and 3D Image Analysis 

Our research is to investigate whether the total number of Mi points mean anything to the object. 
We have discussed Gaussian curvature in the last section. For mean curvatures, Mi also contains 
information about mean curvatures. For instance in a closed digital surface, M4 is independent to 
M3; M4 contains parabolic points (the bended ones). M3 contains the elliptical points since the 
surface point is locally convex. M5 and Mq contain hyperbolic points; the Gaussian curvature is 
negative and the surface point will be locally saddle shaped. 



Lemma 2.3 (a) fcf ' = ^ + V T " f = 4.24913, fc^^^ = ^ - 
kf^ = 0,, k'i"^ = 0; fc("^' = 2^/2 = 2.82843, kf"'' = 0. (c) fcf ^ 

^2^'' t - \/i + f = -0.686875, . (d) kf^ = V¥ = 1.77245, kf^ 
types of Mq digital surface points. 



3.1 Analysis Based on Digital Gaussian Curvatures 



Let 5 be a subset of the 2D manifold M. We define a feature vector called the digital curvature 
vector /s — (7713,7714,7715,7716), where 771^ — \Mi\ with respect to S. We can also split the vector 
into a six component vector to include two cases for M4 and M6 respectively. Basically, the vector 
now contains four components: fs = (l-^^al, IM^I, IM5I, |Mg|). 

The motivation to define such a vector for a local region is based on the corresponding Gaussian- 
Bonnet theorem [T5]. Suppose C is a boundary curve of S that is simply connected and kg is the 
geodesic curvature of C, then 

/ kgdt +11 KdA = 2tt (5) 
Jc J Js 

If C is a 7i-polygon and is the interior angle, then 

(?i-2)-7r+ ( kgdt+ ( ( A'dA^S^^ia, (6) 
Jc J J s 

Therefore, 

9s = 9{fs) = IM3I • A'a + IM5I • ifs + Wh\ ■ K^) ^ J KdA (7) 
can also be used to represent the total geodesic curvature of C : 

[ kgdt ^ 27: ~gs (8) 
Jc 

In image processing, we can define S (or C) as a rectangular region. If S only contains one 
surface point, then gs is just K. In practice, the vector fs can be used for a closed surface or a small 
region centered at one point. Between two regions, there can be intersections or no intersection. 
The region can be 2D or 3D, depending on whether the problem can be projected into 2D easily. 
The most popular shapes of the domain regions are circles/spheres and rectangles/cubes. 

We have used this technique to analyze human facial data. If the size of the region reduced by 
1, 2, 2^, 2'^ times, we will get a sequence of fs, or simple get gs- Then we can see the change 
of the curvatures. Such a method is usually called a multi-scaling method. 

Let us look at the following example. The two original images are shown in Fig 3. 
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Figure 3: Two original human faces from NIST-FRGC data sets 

Fig 4 and Fig 5 show the initial curvature calculation by projection from x— and y— directions. 
)=' — n/2 indicates Ma-points; '+' = —n/2 indicates M5-points; 'X' = — tt indicates Me-points; 



Figure 4: Digital curvature on face-one's surface by projection 




Figure 5: Digital curvature on face-two's surface by projection 

Fig 6. shows the projected Gaussian curvature data for each scale for face-one from 64 x 64 x 64 
to 8 X 8 X 8. Fig 7 is for face-two. In order to be able to observe easily, we will display the same 
size image when the scale changes. 

Analyzing Fig 7 and 8, we saw that the multi-scaling method can identify some of the interesting 
areas of the face. For example, we can identify five areas of interest for face-one and face-two. 
They both contain two areas in the top and in the bottom and one area in the middle. These 
five areas indicate two region on the head-sides, two cheeks, and the nose. The nonzero total of 
Gaussian curvature for a region that remains means that total change has not been canceled in 
this region. 

The result shows that the first person has a flatter face than the second person since the 
method reaches the five interesting areas earlier. We can also identify the flat-or-bending regions 
in the human face as shown in Fig. 9. A total of four such areas can be easily found. 

The advantage of Gaussian curvature-based calculation is that it is not a simple processing of 
pixel averages in the region. The curvature-based calculation is based on geometric and topological 
properties of the 3D object. For instance we know the total Gaussian curvature will be a constant 
as the selected region becomes the whole 3D data array. In summary, the method described in 
this section can identify the five regions of interest and four flat-or-bending regions as the common 
similar characteristics of a human face. In addition, the calculation is much easier than that of 
triangulation based images. The following lemma will provide evidence for this. 

Lemma 3.1 (1) The algorithm designed to obtain scaled local Gaussian curvatures for finding five 
areas of interest on the image of a human face is in 0{nlogn). The space needed is 0{logn). (2) 
There is also an algorithm to solve the same problem in linear space and linear time. 




m o>i i'^) (i) (e) 



Figure 6: Digital curvature scaling for Face-one 

m 0) (c) (d) (e) 



Figure 7: Digital curvature scaling for Face-two 

Proof: (1) Since the 2^ — 1-scaled Gaussian curvature data can be used to calculate the 2'' — 1- 
scaled Gaussian curvature data. We can design a fast linear algorithm due to the reduction of the 
size of the array by half. In such a case, the space needed will also be 0(n) if we just use the input 
array and not the intermediate data. By using the output of the result for each resolution or scale, 
the time complexity will be 0{n\ogn) and the space complexity will be O(logn) (consider only 
the space needed to run the algorithm). (2) If we keep all intermediate data, the space complexity 
will be enlarged to 0{n), but the time complexity will be reduced to 2 • n, where n is the total 
number of points in the array. 

3.2 Investigation Based on Digital Mean Curvatures 

For the mean digital curvature application, we have investigated the calculation based on the 
absolute average for small regions. The geometric meaning of such treatment for digital mean 
curvature is the zigzagged points on the surfaces. Mean curvature zero points indicate the critical 
points that change from inward to outward if it is not a flat point. The result of an example is 
shown in Fig. 9. 

To summarize the findings, Gaussian curvatures describe the global features and the average 
characteristics such as the five regions of a human face, but mean curvatures find local features 
and extreme points such as the nose. 

In the future, we will use this method for analyzing more human facial data. We can normalize 
the image size and calculate the distance for a pair of faces in corresponding scale or nearby 
scales. The following section will present a method for 3D data "rough" classification using digital 
curvature vectors. 

4 3D Object Classification Using Digital Curvature Vectors 

In this section, we will explore a method based on the digital curvature vectors for 3D shape 
similarity analysis and classification of 3D objects. 

There have been many research investigations that use local curvatures in 3D shape similarity 
analysis. Shum, Hebert, and Ikeuchi proposed a method that has 0{n^) time complexity |19| . 




Figure 8: Five areas of interests and four flat-or-bending areas are the common similar charac- 
teristics of human face 
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Figure 9: Digital mean curvature to find a peak point that indicates nose of human face: (a) Image 
is made by a 2 x 2 summation of the absolute values of digital mean curvatures, and (b)Image is 
made by a 4 x 4 summation of the absolute values of digital mean curvatures. 

However, it is not very practical for the purpose of data retrieval. Another disadvantage of the 
continuous curvature method is that the calculation of curvatures return real numbers and that 
introduces precision errors. [TH] [13] [H] A multi-scaling method was proposed to complete this 
task, but it takes more time. 

The technique presented in this section is not to attempt to replace the methods developed 
before. We try to explore more applications using digital curvatures. As we discussed in Section 
2 and 3, the local digital curvature is determined by the local shape of the digital surface. The 
number of each type of surface-point may indicate the features of a 3D object. 

4.1 Feature Vectors of 3D object Based on Digital Curvatures 

First, form a feature vector that only contains the number of digital surface points in each of the 
categories M3, M4, and A/g. 

A 3D digital object may not necessarily be a 3D digital manifold. This is because we have 
strict definitions for 3D manifolds. However, this does not affect the use of curvatures. 

The feature vector is 

fv = ^(iMsl, IM4I, lAfsl, iMel) = (rg, r4, rs, rg) (9) 

where T is the total surface points. 

For example, an object has a total of 1679 digital points. IM3I = 469; IM4I = 995; IM5I = 
183; |Afg| = 32. Because of some non-manifold points, we corrected the data for — M — 's: 484; 



995; 180; 28 ; 240. It also includes the total of 8 none manifold points. Therefore, we get the 
Ratio rs = 0.288267, = 0.592615, = 0.107207, and rg = 0.016677. 
The Euclidean distance of two feature vectors: 

We can also calculate other distances such as the general Minkowski distance. 

d„ = (Et3((^^-2/.r))'/" 

We could also use the scaling method from Section 3 to get more features. Here we use this 
method to do a rough classification for 3D shapes. 

The computing examples are provided in Princeton Benchmark Website, [llj 
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Figure 10: Six 3D objects from Princeton database 



4.2 Similarity and Distance Analysis based on the Feature Vectors 

Feature ratio vectors ex — (r3, r4, r5, r6) where rj. = \Mi^\/T, 
ei = (0.288267,0.592615,0.107207,0.016677) 

62 = (0.262424,0.508752,0.193369,0.044133) 

63 = (0.168149,0.680220,0.144854,0.008895) 

64 = (0.152833,0.711492,0.122506,0.013966) 

65 = (0.148500,0.710425,0.135432,0.007128) 

66 = (0.162700,0.688310,0.140705,0.010093) 
The Distance matrix of two-vector pairs is 



/ 



0.015878586 
0.023580826 
0.032715518 
0.034301844 
0.02609007 





0.041884473 
0.059045308 
0.058376743 
0.04611817 





0.001737666 
0.001390322 
0.000113789 





0.000233753 
0.000980967 



(10) 



0.000727309 / 



In the matrix, describes the distance between object i and object j. The matrix is sym- 
metric so we only present half of it. We can easily see that Objects 1 and Object 2 are closed 
related. Objects 4,5, and 6 are similar. From the data pictures displayed in Fig. 10, we can see 
that is correct. An interesting observation is that Object 3 does not go with three other objects 
in the second category. 
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